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Abstract 

In this paper, in addition to the earlier introduced involutive divisions, we 
consider a new class of divisions induced by admissible monomial orderings. We 
prove that these divisions are noetherian and constructive. Thereby each of them 
allows one to compute an involutive Grobner basis of a polynomial ideal by se- 
quentially examining multiplicative reductions of nonmultiplicative prolongations. 
We study dependence of involutive algorithms on the completion ordering. Based 
on properties of particular involutive divisions two computational optimizations 
are suggested. One of them consists in a special choice of the completion ordering. 
Another optimization is related to recomputing multiplicative and nonmultiplica- 
tive variables in the course of the algorithm. 

1 Introduction 

In paper Q a concept of involutive monomial division was invented which forms the 
foundation of general involutive algorithms 0] for construction of Grobner bases || 
of a special form called involutive. This notion, by a well-known correspondence [[|, [5|] 
between polynomials and linear homogeneous partial differential equations (PDEs) with 
constant coefficients, follows the notion of involutivity for PDEs. An involutive form 
of a system of PDEs is its interreduced completion by the differential consequences 
called prolongations]], incorporating all integrability conditions into the system || ||. 
The integrability conditions play the same role in the completion procedure for PDEs 
as nontrivial S'-polynomials in the Buchberger algorithm |H| |11|] for construction of 
Grobner bases. 

Given a finite polynomial set and an admissible monomial ordering, an involutive 
division satisfying the axiomatic properties proposed in leads to a self-consistent 

1 Prolongation for PDE means its differentiation whereas for a polynomial this means multiplication 
by the corresponding variable. 
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separation of variables for any polynomial in the set into disjoint subsets of so-called 
multiplicative and nonmultiplicative variables. 

The idea of the separation of variables into multiplicative and nonmultiplicative goes 
back to classical papers of Janet || and Thomas [J/j . They used particular separations of 
independent variables for completing systems of partial differential equations to involu- 
tion. More recently one of the possible separations already introduced by Janet || was 
intensively used by Pommaret M in the formal theory of partial differential equations. 
These classical separations allow one to generate the integrability conditions by means 
of multiplicative reductions of nonmultiplicative prolongations. Just this fact was first 
used in as a platform for an involutive algorithm for construction of Pommaret bases 
of polynomial ideals. 

If an involutive division satisfies some extra conditions: noetherity and constructiv- 
ity 0, then an involutive basis may be constructed algorithmically by sequential ex- 
amination of single nonmultiplicative prolongations only. Whereas Thomas and Janet 
divisions satisfy all the extra conditions, Pommaret division, being constructive, is non- 
noetherian. This implies that Pommaret bases of positive dimensional ideals may be 
infinite. The uniqueness properties of involutive bases are investigated in || where a 
special form of an algorithm proposed for construction of a minimal involutive basis 
which is unique much like to a reduced Grobner basis. In addition to the above men- 
tioned classical divisions, in paper [g] two more divisions were introduced which satisfy 
all the extra conditions. 

Recently it was shown [12] that one can also construct different possible separations of 
variables for a fixed monomial set. These separations can not be considered, generally, as 
functions of a set and its element as defined in QTJ] . Nevertheless, the results of paper Jl2[ 



demonstrate for a wide class of divisions how one can change the division dynamically 
in the course of the completion. This increases the flexibility of the involutive technique 
and may also increase the efficiency of computations. 

An involutive basis is a special kind of Grobner basis, though, generally, it may 
be redundant. However, extra elements in the former may facilitate many underlying 
computations. The structure of a Pommaret basis, for example, reveals a number of 



attractive features convenient for solving zero- dimensional polynomial systems [|13| . An 



involutive basis for any division allows one to compute easily the Hilbert function and 
the Hilbert polynomial by explicit and compact formulae [|1^, [L4j] . 

Computation of Janet bases relying upon the original Janet algorithm was imple- 



mented in Reduce and used for finding the size of a Lie symmetry group for PDEs [15 



and for classification of ordinary differential equations admitting nontrivial Lie symme- 



tries |fl6| . The study of algorithmic aspects of the general completion procedure for 
Pommaret division and implementation in Axiom was done in |T7| . The completion 
to involution of polynomial bases for Pommaret division was algorithmized and imple- 
mented in Reduce, first, in [[J, and then with algorithmic improvements in |lj. The main 
improvement is incorporation of an involutive analogue of Buchberger's chain criterion. 



Recently different involutive divisions were implemented also in Mathematica [jM 



In the present paper we introduce a class of involutive divisions induced by admissible 
orderings and prove their noetherity and constructivity. For the new class of divisions, 
along with the classical ones and two divisions of paper we study the stability of 
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the partial involutivity for monomial and polynomial sets under their completion by 
irreducible nonmultiplicative prolongations. We generalize the involutive algorithms to 
different main and completion orderings. The completion ordering serves for selection of 
a nonmultiplicative prolongation to be treated next. In so doing, a completion ordering 
defines the selection strategy in involutive algorithms similar to the selection strategy 
for critical pairs in Buchberger algorithm [|K], p|. For different divisions we find some 



completion orderings which preserve the property of partial involutivity and thereby 
may save computing time. We indicate also a 'pairwise' property which is valid for all 
known divisions. This property can be used to efficiently recompute the separations 
when a new polynomial has to be added. 



2 Background of Involutive Method 

In this section, we recall basic definitions and facts of papers |], which are used in 
the next sections. 

2.1 Preliminaries 

Let N be the set of nonnegative integers, and M = {xf 1 ■ ■ ■ x^™ \ d{ G N} be a set of 
monomials in the polynomial ring R = K[x%, . . . ,x n ] over a field K of characteristic 
zero. 

By deg(u) and degi(u) we denote the total degree of u G M and the degree of 
variable Xi in u, respectively. For the least common multiple of two monomials u, v G M 
we shall use the conventional notation lcm(u,v). If monomial u divides monomial v we 
shall write u\v. In this paper we shall distinguish two admissible monomial orderings: 
main ordering and completion ordering denoted by y and □, respectively. The main 
ordering serves, as usually, for isolation of the leading terms in polynomials whereas the 
completion ordering is used for taking the lowest nonmultiplicative prolongations by the 
normal strategy Jl| and thereby controlling the property of partial involutivity. Besides, 
throughout the paper we shall assume that 

xi y x 2 y ■ ■ ■ y x n . (l) 

The leading monomial and the leading coefficient of the polynomial / Gl with respect 
to y will be denoted by lm(f) and lc(f), respectively. If F C R is a polynomial set, 
then by lm(F) we denote the leading monomial set for F, and Id(F) will denote the 
ideal in R generated by F . The least common multiple of the set {lm(f) \ f G F} will 
be denoted by lcm(F). 

2.2 Involutive Monomial Division 

Definition 2.1 An involutive division L on M is given, if for any finite monomial set 
U C M and for any u G U there is given a sub monoid L(u,U) of M satisfying the 
conditions: 
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(a) If w G L(u, U) and v\w, then v G I/). 

(b) If u, v G 1/ and uL(u, U) H uL(u, 17) 7^ 0, then w G vL(v, U) or v e uL(u, U). 

(c) If v G Z7 and u G 17), then L(v, U) C L(u, 17). 

(d) If y C I/, then L(w, I/) C L(«, V) for all ueK 

Elements of L(u, U), u G U are called multiplicative for w. If w G uL(u, U) we shall write 
and call u (L-)involutive divisor of w. The monomial w is called (L-)involutive 
multiple of w. In such an event the monomial v — w/u is multiplicative for u and the 
equality w = uv will be written as u; = u x v. If m is a conventional divisor of w 
but not an involutive one we shall write, as usual, w — u ■ v. Then v is said to be 
nonmultiplicative for u. 

Definition 2.2 We shall say that an involutive division L is globally defined if for any 
u G M its multiplicative monomials are defined irrespective of the monomial set U 3 u, 
that is, if L(u, U) = L{u). 

Definition |2.1| for every u £ U provides the separation 

{x 1 ,...,x n } = M L (u,U)UNM L (u,U), M L (u,U)nNM L (u,U) = (b (2) 

of the set of variables into two subsets: multiplicative Ml{u, U) C L(u, U) and nonmulti- 
plicative NMl(u, U)nL(u, U) = 0. Conversely, if for any finite set U CM and any u G U 
the separation (0) is given such that the corresponding submonoid L(u,U) of monomials 
in variables in Ml(u, U) satisfies the conditions (b)-(d), then the partition generates an 
involutive division. The conventional monomial division, obviously, satisfies condition 
(b) only in the univariate case. 

In what follows monomial sets are assumed to be finite. 

Definition 2.3 A monomial set U G M is involutively autoreduced or L— autor 'educed if 
the condition uL(u, U) H vL(v, U) — holds for all distinct u,v G U . 

Definition 2.4 Given an involutive division L, a monomial set U is involutive with 
respect to L or L— involutive if 

(Vu G U) (Vw G M) (Bv G U) [uw G f L{v , If) ] . 

Definition 2.5 We shall call the set U„ e [/«M £/ie cone generated by U and denote it 
by C(U). The set U ue u u L(u, U) will be called the involutive cone of 1/ with respect to 
L and denoted by Cl(U). 

Thus, the set U is L— involutive if its cone C(U) coincides with its involutive cone Cl(U). 
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Definition 2.6 An L— involutive monomial set U is called L— completion of a set U C.U 
if 

(Vu G £/) (Vu; G M) (3u G C/) [we IT) ] . 

If there exists a finite L— completion U of a finite set Z7, then the latter is finitely 
generated with respect to L. The involutive division L is noetherian if every finite set Z7 
is finitely generated with respect to L. 

Proposition 2.7 ftjjj If an involutive division L is noetherian, then every monomial 
ideal has a finite involutive basis U. 

Definition 2.8 A monomial set U is called locally involutive with respect to the invo- 
lutive division L if 

(Vu G U) (Vx 4 G NM L (u, U)) (3v G U) [ v\ L (u ■ Xi ) ] . 



Definition 2.9 A division L is called continuous if for any set U G M and for any finite 
sequence {iii}(i<i<k) of elements in U such that 

(Vi < fc) (3xj G NM^^Ui, U)) [ u l+1 \ LUi ■ Xj ] (3) 

the inequality Uj 7^ Uj for z 7^ j holds. 

Theorem 2.10 If an involutive division L is continuous then local involutivity of a 
monomial set U implies its involutivity. 

Definition 2.11 A continuous involutive division L is constructive if for any U C M, 

u G U, Xi G NMl(u, U) such that u ■ Xi has no involutive divisors in U and 

(Vu G £7") (Vxj G NM L (v, U)) (v ■ Xj\u ■ Xi, v • Xj ^ u • Xi) [v ■ Xj G U ue u u L(u, U) } 
the following condition holds: 

(Vtu G U u€U uL(u, U)) [u-Xi^wL(w,UU{w})]. (4) 

Definition 2.12 Let L be an involutive division, and Id(U) be a monomial ideal. Then 
an L— involutive basis U of Id(U) will be called minimal if for any other involutive basis 
V of the same ideal the inclusion U C.V holds. 

Proposition 2.13 f^j IfU CMis a finitely generated set with respect to a constructive 
involutive division, then the monomial ideal Id(U) has a unique minimal involutive basis. 
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2.3 Involutive Polynomial Sets 



Definition 2.14 Given a finite set of polynomials F C E and a mam ordering y, 
multiplicative and nonmultiplicative variables for f E. F are defined in terms of lm(f) 
and the leading monomial set lm(F). 

The concepts of involutive polynomial reduction and involutive normal form are 



introduced similarly to their conventional analogues (Kj [TT(] with the use of involutive 
division instead of the conventional one. 

Definition 2.15 Let L be an involutive division L on M, and let F be a finite set of 
polynomials. Then we shall say: 

(i) . p is L— reducible modulo f G F if p has a term t = au, (a G K\ {0}), u G M such 

that u = lm(f) x v, v G L(lm(f),lm(F)). It yields the L— reduction p —> g = 
p-(a/lc(f))fv. 

(ii) . p is L— reducible modulo F if there is / G F such that p is F— reducible modulo /. 

(iii) . p is in L— normal form modulo F if p is not F— reducible modulo F. 



We denote the L— normal form of p modulo F by NFi(p, F). In contrast, the conven- 
tional normal form will be denoted by NF(p,F). If monomial u is multiplicative to 
lm(f) (/ G F) and h = fu we shall write h = f x u. 

Definition 2.16 A finite polynomial set F is L— autoreduced if the leading monomial 
set lm(F) of F is L— autoreduced and every / G F does not contain monomials which 
are involutively multiple of any element in lm(F). 



Remark 2.17 The further definitions and theorems of this section which involve the 
completion ordering C generalize those in [TJ where C is the same as the main ordering 
y. The proofs of the generalized theorems are immediate extensions of the underlying 
proofs in 0. 



Definition 2.18 An L— autoreduced set F is called (L-) involutive if 

(V/ G F) (Vu G M) [ NF L (fu, F) = ] . 

Given u G M and an L— autoreduced set F, if there exist / G F such that lm(f) C f 
and 

(V/ G F) (Vu G M) (Zm(/) • w C «) [ NFl^/u, F) = ] , (5) 

then F is called partially involutive up to the monomial v with respect to the ordering 
C F is still said to be partially involutive up to v if v C. lm(f) for all / G F. 
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Theorem 2.19 /]!]/ An L—autoreduced set F CM is involutive with respect to a contin- 
uous involutive division L iff the following (local) involutivity conditions hold 

(V/ G F) (Vx, G NM L (lm(f), lm(F))) [ NF L (f ■x i ,F)=0]. 

Correspondingly, partial involutivity (|5p holds iff 

(V/ G F) (Vx, G NM L (lm(f), lm(F))) (lm(f) ■ x, L C v) [ iVF L (/ • x 4 , F) = ] . 



Theorem 2.20 7/F C R is an L— involutive basis ofld(F), then it is also a Grobner 
basis, and the equality of the conventional and L— normal forms NF(p, F) = NFi(p, F) 
holds for any polynomial pGK. If the set F is partially involutive up to the monomial 
v with respect to \Z, then the equality of the normal forms NF(p,F) = NFl(p,F) holds 
for any p such that lm(p) C v. 



Theorem 2.21 Let L be a continuous involutive division, F be a finite L—autoreduced 
polynomial set and NF^(p,F) be an algorithm of L— involutive normal form. Then the 
following are equivalent: 

(i) . F is an L— involutive basis of Id(F). 

(ii) . For all g G F, x G NMi,(lm(g),lm(F)) there is f G F satisfying lm{g) ■ x = 

lm(f) x w and a chain of polynomials in F of the form 

f = /&> /fc-ij • • • > /o> 9o, ■ ■ ■ 1 9m-u 9m = 9 

such that 



NF L (S L (f_ 1 ,f i ),F)=NF L (S(f ,g ),F)=NF L (S L (g j . 1 ,g 1 ),F) = 0. 



where 0<i<k,0<j<m, S(fo,g ) is the conventional S-polynomial / |7Z| [71|/ 
and S^fi^fj) = f ■ x — fj x w is its special form which occurs in involutive 
algorithms j^J. 



Proof (i) =^ (ii) immediately follows from Theorems |2.19| and [2.20| if one takes /o = /, 
go = g. To prove (ii) =^ (i) one suffices to show that NF L (g ■ x,F) = 0. Assume for 
a contradiction that there are nonmultiplicative prolongations which are L— irreducible 
to zero modulo F. Let g ■ x be such a prolongation which is the lowest with respect 
to the main ordering y. This means the partial involutivity of F up to lm(g) ■ x with 
respect to >-. Correspondingly, the condition (ii) implies the representation |TT[ (cf. the 
proof of Theorem 8.1. in IJ) Si,(f : g) = g ■ x — f x w = fi u ij where fi&F and 
Im(fiUij) -< lm(g) ■ x that contradicts NFi(g ■ x, F) ^ 0. □ 
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Corollary 2.22 JJ/ Let F be a finite L— autoreduced polynomial set, and let g ■ x be a 

nonmultiplicative prolongation of g G F . If the following holds 



(V/i G F) (Vu G M) ( lm(h) ■ u C Zm(# • x) ) [ iVF L (/i • u, F) = 

(3/,/o,ff ef) 



lm(f )\lm(f) , lm(g )\lm(g) 
lm(f)\ L lm(g ■ x) , lcm(f ,g ) C ■ x) 
Nmo-m,F)=NF L (g o .M^,F)=0 



i/ien t/ie prolongation g ■ x may be discarded in the course of an involutive algorithm. 



Remark 2.23 Theorem [2.19| is the algorithmic characterization of involutivity whereas 
Theorem [2. 20 relates Grobner bases and involutive bases. Theorem 2.21 and Corol- 



lary |2.22| yield an involutive analog of the Buchberger's chain criterion |10 



Definition 2.24 Given a constructive division L, a finite involutive basis G of ideal 
Id(G) is called minimal if lt(G) is the minimal involutive basis of the monomial ideal 
generated by {lt(f) | / G Id(G)}. 



Theorem 2.25 



A monic minimal involutive basis is unique. 



3 Examples of Involutive Divisions 
3.1 Previously Introduced Divisions 

We give, first, examples of divisions corresponding to separations introduced by Janet, 
Thomas and Pommaret for the purpose of involutivity analysis of PDEs, and two more 
divisions proposed in [0]. For the proof of validity of properties (a)-(d) in Definition |2~I 
for these divisions we refer to [Jj, 0. 

Definition 3.1 Thomas division 0. Given a finite set U C M, the variable xi is 
considered as multiplicative for u G U if degi(u) = max{degi(v) \ v G U}, and nonmul- 
tiplicative, otherwise. 



Definition 3.2 Janet division 0. Let the set U C M be finite. For each 1 < i < n 
divide U into groups labeled by non-negative integers d±, . . . , df 



[di 



{ u G U | dj — degj(u), 1 < j < i }. 



A variable Xj is multiplicative for u G U if i — 1 and degi(u) = max{degi(v) | v G U}, 
or if i > 1, u G [d\, . . . , di-i] and degiiu) = max{degi{v) \ v G [d\, . . . , <ij-i]}. 
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Definition 3.3 Pommaret division ||. For a monomial u = xf 1 --^^ with dk > 
the variables Xj,j > k are considered as multiplicative and the other variables as 
nonmultiplicative. For u = 1 all the variables are multiplicative. 

Definition 3.4 Division I ||. Let U be a finite monomial set. The variable Xi is 
nonmultiplicative for u G U if there is t> 6 U such that 

' ' ' x iZ u = l cm ( u i v ), l <m < [n/2], dj > (1 < j < m) , 
and Xi G {x^, . . . >x^ m }. 

Definition 3.5 Division II [Q. For monomial u = xf 1 • • ■ x d k n the variable Xj is multi- 
plicative if di = d max (u) where d max {u) = max{di, . . . , d n }. 

Remark 3.6 Thomas division, Divisions I and II do not depend on the ordering on 
the variables. Janet and Pommaret divisions, as defined, are based on the ordering 
given in ([!]) . Pommaret division and Division II are globally defined in accordance with 
Definition |2.1| 

All these divisions are constructive, and except Pommaret division they are noetherian 

EE- 

3.2 Induced Division 

Now we consider a new class of involutive divisions induced by admissible monomial 
orderings (cf. [0). 

Definition 3.7 Induced division. Given an admissible monomial ordering >-0 a variable 
Xi is nonmultiplicative for u G U if there is v G U such that v -< u and degi(u) < degi(v). 



Proposition 3.8 The separation given in Definition \3. Tj is an involutive division. 

Proof Let L y (u, U) be the submonoid generated by multiplicative variables. We must 
prove the properties (b-c) in Definition ^71] because (a) and (d) hold obviously. 

(b) Let there be a monomial w such that w G uL^(u, U) fl vL y (v, U) with u,v G U 
and u^v. Assume u >- v and —<v\u. Then, there is a variable x\(lcm(u,v)/u) such that 
x ^ L^(u, U). Since v\w we obtain x\(w/u) that contradicts w G uL y (u, U). Thus, v\u 
and w = v x (w/v) = v x [{w /u){u/v)}. This yields u G vL y {y, U). 

(c) Let v G uL^{u,U) and w G vL y (v,U) with u,v G U, and, hence, u\v and u\w. 
Suppose w G" uL y (u, U). It follows the existence of a variable x\(w/u), ->x\{v/u) and a 
monomial t~<u^v,tEU such that x\(lcm(u, t)/u). This suggests that x\(lcm(v, t)/v) 
at i -< v , contradicting our initial assumption. □ 

2 This ordering is generally different from the main ordering introduced in Sect. 2.1. 
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Remark 3.9 Generally, an ordering y defining Induced division implies some variable 
ordering which is not compatible with (jl]). However, below we assume that the ordering 
y is compatible with ([I]). 

To distinguish the above divisions, the abbreviations T, J, P, I, II, will be sometimes 
used. For illustrative purposes we consider three particular orderings to induce involutive 
divisions: lexicographical, degree-lexicographical and degree-reverse- lexicographical. To 
distinguish these three orderings we shall use the subscripts L, DL, DRL, respectively. 
There are certain relations between separations generated by those divisions. 

Proposition 3.10 For any U , u G U and y the inclusions Mt{u,U) C Mj(u,U), 
Mt(u, U) C Mi(u, U), Mt(u, U) C Mn y {u, U) hold. IfU is autoreduced with respect to 
Pommaret division, then also Mp{u,U) C Mj(u,U). 

Proof The inclusion Mt{u,U) C Md^(u,U) follows from the observation that x G 
T(u, U) implies x G D y {u, U). The other inclusions proved in Jl], 0. □ 

The following example explicitly shows that all eight divisions we use in this paper are 
different. In the table we list the multiplicative variables for every division. 



Example 3.11 (|TJJ] Multiplicative variables for elements in the monomial set U = 

{x 2 y,xz,y 2 ,yz, z 3 } (x y y y z) for different divisions: 



Monomial 


Multiplicative variables 


T 


J 


p 


I 


II 


D L 


D DL 


Ddrl 


2 

x y 


X 


x,y,z 


y,z 


X 


X 


X 


X 


X 


xz 




y,z 


z 


X 




X 






y 2 


y 


y,z 


y,z 


y 


y 


x,y 


x,y 


y 


yz 




z 


z 




y,z 


x,y 


x,y,z 


x,y,z 


z 3 


z 


z 


z 


z 


z 


x,y,z 


z 


z 



Proposition 3.12 Induced division is noetherian, continuous and constructive. 

Proof Noetherity. follows immediately from noetherity of Thomas division and the 
underlying inclusion in Proposition |3.10| . 

Continuity. Let U be a finite set, and {«j}(i<j<Af) be a sequence of elements in U 



satisfying the conditions (y). In accordance with Definition [2.9| we shall show that there 
are no coinciding elements in the sequence for each of the two divisions. There are the 
following two alternatives: 

(i) Ui = Ui-i ■ Xj\ (ii) ui ^ Ui-i ■ xj . (6) 
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Extract from the sequence {ui} the subsequence {t k = u ik \ (i<k<K<M) of those elements 
which occur in the left-hand side of relation (ii) in (^|). 

Show that t k+1 \ D lcm(t k+1 ,t k ) and ~<t k \t k+1 . We have t k+1 x w k+1 = u ik ^ ■ x jk = 
t k ■ v k where ^w k+ i\v k . Indeed, suppose w k+ i\v k . Apparently, we obtain the relation 
t k+ \ — m ■ zi where % k < I < i k +i, and the variable Xj t G NMd(ui,U), which figures 



in Definition |2]9| of the sequence satisfies Xj^Wk+i and -iXj t \zi. This suggests, by 

definition of the division, the existence ofp G U such that p -< ui and deg^ (ui) < deg^ (p) . 
Since p -< t k+1 and deg^ui) = deg^ , it contradicts multiplicativity of Xj t for t k+1 . 
Therefore, we obtain the relation 



| t k -v k = t k+ i x w k+1 , 

\ gcd(v k} w k+ i) = gcd(v k ,w k ) = 1, 

which, by Definition |3.7| , implies t k y t k+ i since w k ^ 1 for all fc. 

It remains to prove that elements in the sequence {wj}(i<j<M) which occur in the left- 
hand side of relation (i) in (j^) are also distinct. Assume for a contradiction that there 
are two elements Uj = u k with j < k. In between these elements there is, obviously, 
an element from the left-hand side of relation (ii) in (|B|). Let w im (j < i m < k) be 
the nearest such element to Uj. Considering the same nonmultiplicative prolongations 
of u k as those of Uj in the initial sequence, one can construct a sequence such that 
the subsequence of the left-hand sides of relation (ii) in @ has two identical elements 

Constructivity. Let u,Ui G U, v G D y (ui,U) and Xj G" D y (u,U) be such that 
u ■ Xi = U\V x w, w G D y (uiv,U U {wit>}). Show that w G D y (ui,U). Assume that 
there is Xj G" D y (u\,U) satisfying Xj\w. This implies the existence t G f/ satisfying 
£ -< w, degj(t) > degj(u\). Then, because t -< U\V, the condition — ix^ | w leads to the 
contradictory condition degj(t) > degj(uiv). Therefore, Xj\v. □ 



4 Completion of Monomial Sets to Involution 

If U is a finitely generated monomial set with respect to the fixed involutive division L, 
then its finite completion gives an involutive basis of the monomial ideal generated by 
U. There may be different involutively autoreduced bases of the same monomial ideal. 



For instance, from Definitions |3.1| and |3.2| it is easy to see that any finite monomial set 
is Thomas and Janet autoreduced. Therefore, enlarging a Thomas or a Janet basis by 
a prolongation of any its element and then completing the enlarged set leads to another 
Thomas and Janet basis, respectively. Similarly, Division I and Induced division do not 
provide uniqueness of involutively autoreduced bases whereas Pommaret division and 
Division II do, as well as any globally defined division ||. 

4.1 Completion Algorithm 

Theorem 4.1 If U is a finitely generated set with respect to a constructive involutive 
division, then the following algorithm computes the uniquely defined minimal completion 
U of U , that is, for any other completion U the inclusion U C U holds. 
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Algorithm InvolutiveCompletion: 

Input: U, a finite monomial set 
Output: U, an involutive completion of U 

begin 1 

U :=U 2 

while exist it G U and x G NMl(u, U) such that 3 

u ■ x has no involutive divisors in U do 4 

choose any C and such it and x with the lowest it ■ x w.r.t. C 5 

U:=Uu{u-x} 6 

end 7 

end 8 



Proof This completion algorithm is a slightly generalized version of that in [I] where 
ordering C is assumed to be fixed in the course of the completion. As proved in |l] (see 
the proof of Theorem 4.14), the output U of the algorithm and the number of irreducible 
prolongations are invariant on the choice of ordering in line 5. 

Corollary 4.2 If set U is conventionally autoreduced, then the algorithm computes the 
minimal involutive basis of monomial ideal Id(U). 

In practice, in the course of the completion one has to choose the lowest nonmultiplicative 
prolongation and to check whether it has an involutive divisor in the set. The timing of 
computation is thereby determined by the total number of prolongations checked. 

Theorem 4.3 The number of nonmultiplicative prolongations checked in the course of 
algorithm InvolutiveCompletion with a constructive division L is invariant on the 
choice of completion ordering in line 5. 



Proof As it has been noticed in the proof of Theorem |47l|, the number of irreducible 
prolongations, as well as the completed set itself, is invariant. Therefore, we must prove 
invariance of the number of reducible prolongations. 

Let there be two different completion procedures of U to U based on different choice 
of completion orderings. Assume that the first procedure needs more reducible pro- 
longations to check than the second one. Let u ■ x = v x w (u,v G U\) be the first 
prolongation checked in the course of the first procedure and such that in the course of 
the second one the prolongation is not checked. This suggests x G Ml(u, U?) where U2 
is the current set for the second procedure. If w ^ 1, by admissibility of a completion 
ordering, we obtain u x x — v x w (u, « G U2). From property (b) in Definition [2.1| we 
deduce u G vL(v,U<i), and, hence, x cannot be nonmultiplicative for it as we assumed 
for the first procedure. 

If w = 1 we find that x G NMl(u, Ui)C]Ml(u, JJq) where u and v = ux are elements in 
both Ui, U 2 - From the property (d) in Definition |2 . 1| and invariance of the final completed 
set U is follows that in some step of the second procedure x becomes nonmultiplicative 
for it. Then the prolongation it • x will be also checked that contradicts our assumption. 

□ 
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This theorem generalizes Remark 3.13 in [|14[] which is concerned with L— autoreduced 
sets and fixed completion orderings. 

Example 4.4 (Continuation of Example |3.11| ). The minimal involutive bases of the 
ideal generated by the set U = {x 2 y, xz, y 2 , yz, z 3 } (x y y y z) are given by 



TT r223 222 22 2223222 2 23 

Ut = {x y z ,x y z ,x y z,x y ,x yz ,x yz ,x yz,x y,x z , 

111 2 3 2 2 2 2 3 2 32 

x z ,x z,xy z ,xy z ,xy z,xy , xyz , xyz , xyz, xz ,xz , 

2 3 222 2 3 2 3i 

xz,y z ,y z ,y z,y ,yz ,yz ,yz,z \ , 
Uj = {x 2 y,x 2 z,xy 2 ,xyz,xz,y 2 ,yz,z 3 } , 
Up = {x 2 y,x 2 z,xy 2 ,xyz,xz,y 2 ,yz,z 3 , . . . ,x k y, . . . ,x l z, . . .} , 

tt r22322222 2223222 2 23 

(J i = {xyz,xyz,xyz,xy,x yz ,x yz ,x yz,x y,xy z , 
xy 2 z 2 , xy 2 z, xy 2 , xyz 3 , xyz 2 , xyz, xz 3 , xz 2 , xz, y 2 z 3 , y 2 z 2 , 
y 2 z,y 2 ,yz 3 ,yz 2 ,yz,z 3 } , 

f~T r222 2 2 3l 

U n = {x y , x y, xy , xyz, xz, y ,yz,z }, 

Ft r2 2 22 3t 

U L = {x y,xz ,xz,y ,yz ,yz,z } , 
U DL = {x 2 y,xz,y 2 ,yz,z 3 } , 
U DRL = {x 2 y,xy 2 ,xz,y 2 ,yz,z 3 } , 

where k, I G N (k, I > 2), and subscripts in the left-hand sides stand for different invo- 
lutive divisions considered in Section 3. This example explicitly shows that Pommaret 
division is not noetherian, since it leads to an infinite monomial basis. 



4.2 Pair Property 

In the course of algorithm InvolutiveCompletion the current monomial set U is en- 
larged by irreducible nonmultiplicative prolongations in line 6. As this takes place, for 
a nonglobally defined division one has to recompute the separation into multiplicative 
and nonmultiplicative variables for all monomials. The next definition and proposition 
give a prescription for efficient recomputing. 

Definition 4.5 We shall say that an involutive division L is pairwise if for any finite 
set U and any u E U (U \ {u} ^ 0), the following holds: 

L(u, U) = n veU \ {u} L(u, {v}) 

or, equivalent ly, 

M L (u, U) = n veU \ {u} M L {u, {«}), NM L (u, U) = U veU \ {u} NM L (u, {«}). (7) 



Therefore, for a pairwise division L and a monomial set U the correction of the separation 
due to enlargement of U by an element v is performed by formula 

NM L (u, U U {v}) = NM L (u, U) U NM L (u, {u, v}). (8) 
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Proposition 4.6 All the above defined divisions are pairwise. 



Proof Pommaret division and Division II, as globally denned divisions, are trivially 
pairwise. 

Thomas division. Since 

max{degi{w) \ w G U} = max v( zu\{ u }{max{degi(u),degi(v)}}, 



Definition [3J] implies apparently (|7])- 

Janet division. For i — 1, by our convention (H) and Definition 3.2 , Janet case is re- 



duced to Thomas one, and we are done. Let now i > 1 and deg\{u) = d±, . . . , degi~i(u) = 
di-i. If the group [d\, . . . , di-i] of elements in U contains, in addition to u, some extra 
elements, then 

max{degi(w) \ w G [d x , . . . , cfi-i]} = max ve [ dlr .. tdi _ 1 ]\{ u }{max{deg i (u),deg i {v)}}, 

and Xi G Mj(u, U), otherwise. This suggests the pairwise property. 

Division I and Induced division. For these division the pairwise property follows 



immediately from Definitions [3.4| and |37?|- □ 



4.3 Monotonicity 

Consider now another optimization related to the choice of a nonmultiplicative prolon- 
gation in line 5 of algorithm InvolutiveCompletion. The choice of the lowest prolon- 
gation with respect to some fixed ordering C is called normal selection strategy 

Definition 4.7 Given a division L and an admissible ordering C, a monomial set U 
will be called complete up to monomial w with respect to C if 

(Vw G £/) (Vx G NM L (u, U))(u-xQw) [u-xe C L (U) } , (9) 



where Cl(U) is involutive cone of U by Definition |2.5| . We call monomial w bound of 
completeness for f/. If u • x □ w for all u G U, x G NM L (u, U), then we shall still say 
that U is completed up to to. 



Definition 4.8 We shall call division L monotone for C if for any set U and any 
monomial w G M satisfying (P) the following holds: 

(Vu G C/) (Vx G NM L (v, U)) (vx& C L (U)) [U\J{v ■ x} is complete up to w }. 

We shall say that L is monotone if its monotonicity holds for any ordering C 

Thus, monotonicity means that enlargement of U by an irreducible nonmultiplicative 
prolongation does not decrease its completeness bound. 
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Remark 4.9 If a division L is monotone for an ordering C, then the choice of the latter 
as a completion ordering is beneficial for the algorithm InvolutiveCompletion. By 
Theorem ( |4.3| ), the total number of prolongations checked is invariant on the ordering. 
Monotonicity of the latter allows one to omit recomputing separations and checking 
prolongations which are lower than the current completeness bound. 

Now we consider the monotonicity properties of different divisions defined in Sect. 2. 
Pommaret division and Division II, as globally defined, are trivially monotone. 

Proposition 4.10 Thomas division is monotone. 

Proof From Definition |3.1| it follows immediately that T(u, U) = T(u, U U {v ■ x}) for 

any v G U, x G NM T (v, U). □ 

Proposition 4.11 Janet division is monotone for lexicographical ordering. 

Proof Denote the lexicographical completion ordering compatible with ([I]) by \Zlbx- 
Consider a nonmultiplicative prolongation v-Xj G" Cj(U) (v G U) such that v-Xj Z\iex w 
where w is the completeness bound of U in accordance with (|S]). 
Suppose there is a pair {u G U,Xk}, satisfying 

UX k eCj(U), U ■ X k Cj(UU {vXj}), U ■ X k ^Lex W \Z L ex V ■ Xj, (10) 

and consider the lowest such pair with respect to \ZLex- 
If Xk G J(u, U) we obtain 

x (u) = degi(v) + 1, deg k {u) < deg k {v) if j = 1, 

i(u) = degi(v) (i < j), deg^u) = degj(v) + 1, deg k {u) < deg k (v) if j > 1. 

Here k > j and, if k — j > 1, then deg m (u) = deg m (v) for alH < m < k. Consider now 
two alternatives: 

(i) w G U. In this case conditions fllOD are contradictory since from the rightmost 
condition it follows degiiu) = degiiw) (i < k) and deg k {w) > degk(u), that is, Xk G 
NMj{u,U). 

(ii) w G" U . Then there is t G U such that w G tJ(t, U). Because Xk G J(u, U), for 
some 1 < p < k we have degiit) = degiiu) = degiiyxj) where i < p and deg p+ i(t) < 
deg p+ i(u) < deg p+ \(w) . Thus we obtain contradiction with x p+ \ G J(t, U) which follows 
from w G tJ(t, U). 

It is remains to prove that if Xk G NMj(u, U), then u ■ Xk G Cj(U U {v ■ Xj}). If 
u ■ x k G U we are done. Otherwise, we have u ■ Xk = <?i ?"i for some q x G U, T\ G J(q, U) 
and 7*1 G" J(C/ U {t> • x^}). Hence, there is x^ri, x ix G NM J (q 1 ,U U {t> • x^}), and 
deg(qi ■ x^) \Ziex deg(u ■ x k ). Then, by our assumption that prolongation u ■ Xk is the 
lowest satisfying (|T0|) , we have q\ ■ x^ = q2 X r 2 , g2 G C/, r 2 G J(g2, ?7 U u • xj). By 



property (d) in Definition pTT] , it yields r 2 G J(?2, C^), and, hence, q\ x = g 2 X r 2 in 



{7. This is impossible, because any monomial set is Janet autoreduced. □ 

Remark 4.12 Janet division is not monotone for degree-lexicographical and degree- 
reverse-lexicographical orderings as the following example shows. 
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Example 4.13 Consider the conventionally autoreduced set U = {xz 2 ,x 2 z,yzt 2 }. Let 
completion ordering Z be degree-lexicographical or degree-reverse-lexicographical or- 
dering with x Z y Z z □ 1 £7 is complete up to w = x 2 ?/^. The lowest irreducible 
prolongation is xyzt 2 □ to. The next one in the set {xz 2 , x 2 z, yzt 2 , xyzt 2 } is xyz 2 IZ w. 



Example 4.14 Consider the set {7 = {xy 2 u; 2 , xzt, yzt} and Division I generating the 
separation: 



Monomial 


Division I 


M I 


NMi 


p = xy 2 w 2 


x,y,w 


z, t 


v = xzt 




y,w 


u = yzt 


y,z,t,w 


X 



Let xZlyZlzZltZlw and □ be any of the orderings: lexicographical, degree- 
lexicographical or degree-reverse-lexicographical. The bound of completeness for U is 
xyw 2 . We find that v ■ w = u ■ x = xztw is the lowest irreducible prolongation in U, and 
the next one for U U {xztw} is u ■ w Z xyw 2 . Therefore, Division I is not monotone for 
three orderings considered. 



Proposition 4.15 Induced division is monotone for the ordering which induces this 
division. 

Proof By Definition |3.7| of D n , enlargement of U by irreducible nonmultiplicative 
prolongation v • Xj Z w (v 6 U) does not change the reducibility properties of those 
prolongations uxk (u G U) which satisfy uxu Z v ■ Xj. □ 



5 Construction of Involutive Bases for Polynomial 
Ideals 

In this section we present the following algorithm for computation of minimal involutive 
bases of polynomial ideals which generalizes the algorithm of paper to different 
completion and main orderings. 

Algorithm Minimallnvolut i veBasis : 

Input: F, a finite polynomial set; L, an involutive division; 

y, a main ordering; Z, a completion ordering 
Output: G, the minimal involutive basis of Id(F) if algorithm terminates 
begin 1 
F := Autoreduce(F) 2 
choose g G F with the lowest lm(g) w.r.t. -< 3 
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T:={{g,lrn(g),(t>)}; Q := 0; G := {#} 4 

for each f E F \ {g} do 5 

Q:=gu{(/,M/)>0)} 6 

repeat 7 

ft, := 8 

while Q ^ and ft = do 9 

choose g in (g,u,P) G Q with the lowest lm(g) w.r.t. -< 10 

g : =g\{((7,n,p)} " n 

if Criterion(g,u,T) is false then h:=NFi(g,G) 12 

if ft^O then G:=GU {ft} 13 

if lm{h) = lm(g) then T := T U {(ft, u, P n NM L (h, G))} 14 

else T:=TU {(ft, Zm(ft), 0)} 15 

for each (/, u, £>) G T s.t. lm(f) y lm(h) do 16 

T:=T\{(f,v,D)}- Q:=QU{(f,v,D)}; G:=G\{f} 17 

for each (/, v,D) G T do 18 

T := T \ {(/, v, D)} U {(/, v,D(~) NM L (f, G))} 19 

while exist (g, u, P) G T and a; G NM L (g, G)\P and, if Q ^ 0, 20 

s.t. Zm(p • x) -< lm{f) for all / in (/, v,D) G Q do 21 

choose such (g,u,P),x with the lowest lm(g) ■ x w.r.t. C 22 

T:=T\{(g,u,P)}U{(g,u } PU{x})} 23 

if Criterion(g ■ x, u, T) is false then ft := NF^g ■ x,G) 24 

if ft^O then G:=GU {ft} 25 

if Zm(ft) = lm(g ■ x) then T := T U {(ft, u, 0)} 26 

else T := TU {(ft, Zm(ft), 0)} 27 

for each / in (/, v,D) G T with Zm(/) ^ Zm(ft) do 28 
T:=r\ {(/,«,£))}; Q := Q U {(/, t;, £>}); G :—G \ {/} 29 

for each (/, u, D) G T do 30 

T :—T \ {(/, f , £))} U {(/, «,Dfl iVM L (/, G))} 31 

until g^0 32 

end 33 

Criterion(g, u, T) is true provided that if there is (/, v,D) £ T such that lm(f)\Llm(g) 
and lcm(u, v) C lm(g). Correctness of this criterion, which is just the involutive form Q 
of Buchberger's chain criterion [ TO | , is provided by Corollary |2.22| . 



Theorem 5.1 Let F be a finite subset ofR and L be a constructive involutive division. 
Suppose the main ordering y is degree compatible. Then the algorithm Minimallnvo- 
lutiveBasis computes a minimal involutive basis of Id(F) if this basis is finite. If L is 
noetherian, then the basis is computed for any main ordering. 



Proof The proof is the same as in || and based on Theorems |2.19| , [2.20| and 4.1 



Corollaries ^2| and gj. □ 

Proposition 5.2 The conventional autoreduction of the input polynomial set in line 2 
is optional and may be omitted. 
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Proof Let F be a non-autoreduced set and the algorithm start with line 3. Subsequent 
to the initialization in lines 4-6 the upper while-loop selects, first of all, those polyno- 
mials in the triple set Q which have the same leading term as the element in G = {g}. 
If there is such a polynomial in the triple set Q with nonzero involutive normal form 
h computed in line 12, then lm(h) -< lm(g). It follows from lines 13 and 17 that G 
becomes the one-element set {h} as an input for the lower while-loop. 

Thus, by restriction in line 21 for nonmultiplicative prolongations checked and redis- 
tribution of polynomials in line 29, in every step of the algorithm we have lm(g) -< lm(f) 
for any g in (g, u, P) G T and / in (f,v,D) G Q whenever the set Q is nonempty. 

Furthermore, as proved in |], 0, in some step of the algorithm a polynomial h is 
added to the current polynomial set G in line 13 or in line 25, such that h is an element 
in the reduced Grobner basis of Id(F) with the lowest leading monomial with respect 
to the main ordering -<. It implies the reduction of G to the one-element set G = {h}, 
and transfer of the rest to Q. Then G is sequentially completed by other polynomials 
from the reduced Grobner basis and their nonmultiplicative prolongations. In so doing, 
the completion of lm(G) due to the redistribution of polynomials between sets T and Q 
in lines 17 and 29 is monotone with respect to -<;. 

Therefore, the output of algorithm MinimallnvolutiveBasis irrespective of au- 
toreduction in line 2 is the same as it would be for the reduced Grobner basis in the 
input. □ 

Remark 5.3 The choice of a completion ordering which is monotone for L preserves, 
obviously, the partial involutivity of the intermediate polynomial set G in the case of 
its enlargement in line 25, if lm(h) = lm(g ■ x). Therefore, similar to the monomial case 
(c.f Remark |4.9| ), this saves computing time for recomputing separations and checking 
irreducibility of nonmultiplicative prolongations unless L is globally defined anyway. 



6 Conclusion 

The above described optimizations concern only that part of computing involutive bases 
which is related to completion by nonmultiplicative prolongations with irreducible lead- 
ing terms. Another important step is to search for an involutive divisor among the 
leading monomials of an intermediate basis. This is important for efficient computation 
of the involutive normal form in lines 12 and 24 of algorithm MinimallnvolutiveBa- 



sis. Some related optimizations are considered in |rjj for the purpose of implementing 
the algorithm InvolutiveCompletion in Mathematica for divisions of Sect. 3. 

A promising way to the further optimization of computation is related to the ideas of 



paper |12| . By appropriate dynamical refinement of an involutive division in the course 
of computation, one can decrease the total number of nonmultiplicative prolongations 
to be checked. This may lead to a notable reduction of computing time. 

Algorithm MinimallnvolutiveBasis has been implemented in Reduce for Pom- 
maret division. Computer experiments showed that this algorithm is somewhat faster 
than our previous version of involutive algorithm also implemented in Reduce for Pom- 
maret bases ffl. For a nonglobally defined division the difference in speed is to be much 
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greater as algorithm MinimallnvolutiveBasis deals with fewer intermediate polyno- 
mials and avoids intermediate autoreductions [|f. 

With the new implementation one needs, for example, 57 seconds to compute a 
degree-reverse-lexicographical Pommaret basis for 6th cyclic roots on an Pentium 100 
Mhz computer, and 30 seconds for the 6th Katsura system. By comparison, the PoSSo 
software for computing Grobner basesQ needs for these examples 24 and 36 seconds, 
respectively. 
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